% PURPOSE: An example using barmap
%          to examine the scatter of a
%          variable using a map
%---------------------------------------------------
% USAGE: barmap_d
%---------------------------------------------------

clear all;
% A spatial dataset on crime, household income and housing values
% in 49 Columbus, Ohio neighborhoods
% from:
% Anselin, L. 1988. Spatial Econometrics: Methods and Models,
% (Dorddrecht: Kluwer Academic Publishers).

% 5 columns:
% column1 = crime
% column2 = household income
% column3 = house values
% column4 = latitude coordinate
% column5 = longitude coordinate
load anselin2.data;
crime = anselin2(:,1);
latt = anselin2(:,4);
long = anselin2(:,5);
n = length(latt);
crime = roundoff(crime,0);

load columbus.poly; % see the file for polygon format
carte = columbus;

barmap(latt,long,crime,carte,[],1);